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We consider surface diffusion of a single particle, which performs site-to-site under-barrier hopping, 
fulfils intrasite motion between the ground and the first excited states within a quantum well, 
and interacts with surface phonons. On the basis of quantum kinetic equations for one-particle 
distribution functions obtained earlier we study the coherent and incoherent motion of the adparticle. 
In the latter case we derive the generalized diffusion coefficients and study various dynamic regimes 
of the adparticle. The critical values of the coupling constant Gcr(r, which separate domains 
with possible recrossing from those with the monotonic motion of the adparticle, are calculated as 
functions of a temperature T and a vibrational frequency Q.. These domains are found to coincide 
with the regions where the experimentally observed diffusion coefficients change its behavior from 
weakly dependent on T to quite a sensitive function of the temperature. We also evaluate the off- 
diagonal (relative to the site labels) distribution functions both in the Markovian limit and when the 
memory effects become important. The obtained results are discussed in the context of the "long 
tails" problem of the generalized diffusion coefficients, the recrossing/multiple crossing phenomena, 
and an eventual interrelation between the adparticle dynamics at short times and the temperature 
dependence of the difi'usion coefficients measured experimentally. 
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I. INTRODUCTION 



A phenomenon of quantum diffusion of light particles 
(mostly, hydrogen and its isotopes), adsorbed at solid 
surfaces, in the last decades is a subject of joint efforts 
of investigators in various fields of sciences. An attention 
of the scientists is dictated by its relevance in the tech- 
nological processes like heterogeneous catalysis fuel 
cell production [2|], chemical reactions of hydrogen trans- 
fer Q , and series of physical phenomena occurring at the 
fluid-gas interfaces. On the other hand, the diffusion of 
hydrogen is of fundamental interest from a basic physi- 
cal point of view, being a favorite system for theoretical 
analysis and computer simulations [3, B ■ Wide per- 
spectives have been opened after creation of the scanning 
tunnelling microscope (STM) Q, allowing to "touch" a 
single molecule at the surface and to carry out experi- 
ments which previously were only imagined. At the same 
time a rapid development of powerful methods of com- 
puter experiments like quantum molecular dynamics (loj 
or Monte Carlo Wave Function formalism ll| allows a 
direct analysis of the particle trajectories in real space 
and time. Recent results showing that even such the 
"heavy" atoms like Cr on the Au(lll) surface [11] or 
Na on Cu(OOl) surfaces [Tsj manifest a great deal of the 
underbarrier tunnelling bring us to a conclusion that a 
traditional viewpoint on the quantum diffusion as an in- 
herent feature of light particles only is far from reality, 
and a fresh look at such processes is quite topical. 

A task for the theorists - to create reliable models de- 
scribing a quantum surface diffusion and to develop ef- 
fective methods of the calculation of diffusion coefficients 
with taking into account all the interactions ( "adsorbate- 
adsorbate", "substrate-adsorbate" etc.) has been ac- 



complished to a great extent. 
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tioned the works '14l4l6l| in which a bulk quantum dif- 
fusion of light particles has been studied: its descrip- 
tion requires similar theoretical methods, and elaborated 
schemes can be considered as starting points for investi- 
gation of the surface quantum diffusion. They were fol- 
lowed by the papers where quantum hopping was moved 
from the bulk to the surface 
concept of a small polaron 
tion [6] has been applied with going beyond the linear 
"adatom-phonon" coupling to consider the anharmonic 
terms in "adsorbate-substrate" interaction [l^ . The lat- 
ter case together with consideration of electronic friction 
in the system [ij, and direct "adsorbate-adsorbate" 
interaction [2l| are very important, because they provide 
additional channels of particle scattering and ensure fi- 
nite values of the diffusion coefficients 22] . 

The diffusion coefficients are traditionally determined 
via the Green-Kubo relations ^22] , the low- frequency and 
small wave- vector limit of the dynamic structure factor 
[13 J or zeroth moments of the "flux-flux" time corre- 
lation functions for classical [l^, [13] or quantum systems 
[3, [B]- The latter approach is of particular interest be- 
cause it allows one not only to study in detail the well- 
known "recrossing/miltiple crossing" problem [3, [1], but 
also to derive a new quantum transition state theory [25| . 

The quantum diffusion coefficient is known to consist 
of two terms of a different physical origin. A coherent 
term [l^ characterizes the way in which the dephas- 
ing limits the band motion of the adatom by destruction 
of the coherence of the hopping probabilities when the 
adatom-thermal bath coupling induces random fluctu- 
ations of each phase. This term has a pure quantum 
nature and is related to the competition between the 
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tunnelling mechanism, which tends to preserve the co- 
herence, and the dephasing which characterizes damping 
due to the scattering process. The coherent part Dcoh{T) 
of the diffusion coefficient weakly depends on tempera- 
ture at low T in contrast to the incoherent one Din{T), 
which decreases to zero when T — > 0. The incoherent 
contribution describes processes in which the surface dy- 
namics induces fluctuations of the tunnelling matrix el- 
ements between two Wannier states allowing the adpar- 
ticle to perform a transition from one state to another 
by creating or annihilating surface phonons. Though the 
temperature dependence of the diffusion coefficients was 
studied profoundly for both small polaron model jl8j and 
its modification [6|, some questions remain unclarified 
(especially, in the weak-coupling limit, where the contri- 
butions of Dcoh{T) and Din{T) are of the same order). 
The most challenging are the following: i) a justification 
of the multiphonon expansion; ii) a correct definition of 
the activation energy (e.g. when the main contribution 
arises from the acoustic branch of phonon spectrum) ; iii) 
an introduction of the additional channels of particles 
scattering ensuring correct values of the diffusion coef- 
ficients; iv) an influence of the memory effects on the 
temperature dependence of the diffusion coefficients and 
on the short-time dynamics of the adsorbate. 

The latter problem should be considered from several 
standpoints. First of all, a thorough study of the short- 
time dynamics of the adparticle allows one to distin- 
guish between various scenarios of the adsorbate motion 
(a transition from the coherent to incoherent regime, a 
presence of multiple or long hopping etc.) that could be 
helpful for a deeper insight into the microscopic picture 
of the process [26]. Secondly, we can answer the ques- 
tion: what can enhance or suppress the particle motion. 
For instance, an eventual recrossing reduces the values 
of the diffusion coefficients, whereas a multiple crossing 
increases them. At last, such a theoretical analysis can 
give some recommendations for experimentalists how to 
perform evaluation of the diffusion coefficients more ef- 
fectively. It is known 9] that at temperatures above 80 K 
the diffusion rate of hydrogen is too fast to be followed by 
the standard atom-tracking technique, while below 50 K 
the opposite problem occurs, and it is necessary to min- 
imize the influence of the STM tip on the adsorbate due 
to a prolonged interaction. Thus, it would be tempting to 
relate a change in the character of the short-time dynam- 
ics of the adsorbate to possible crossover from one typical 
temperature behavior of the diffusion coefficients to an- 
other and to give a prognosis about the T-dependence of 
the diffusion coefficient having only an information about 
the adatom dynamics at the initial stage of its motion. 

In the present paper, which is a logical continuation 
of Ref. [26|, we try to give an answer to the question 
about the interrelation between the process of the deco- 
herence in the "adsorbate-substrate" system (leading to 
the dissipative dynamics of the adsorbate) and the mem- 
ory effects. Here we also use the method of the quantum 
kinetic equations p^ [22l.[27t. 



The subject of our study is a single adsorbate which 
performs an underbarrier hopping to the nearest adsorp- 
tion sites, fulfils intrasite motion between two different 
quantum states within a quantum well, and interacts 
with acoustic phonons. We define the conditions under 
which the adatom dynamics is definitely nonlocal in time, 
derive expressions for the generalized (time-dependent) 
diffusion coefficients, and study the influence of their 
long-time asymptotics on the temperature behavior of 
the transport coefficients measured experimentally. A 
particular emphasis is put on the critical regimes sep- 
arating dynamics with mainly coherent contribution to 
the diffusion coefficient from the dominance of that of 
the incoherent origin. We show that a transition domain 
from the oscillating dynamics of the adparticle to the 
monotonic motion subject to the critical coupling con- 
stant Gcr{T,fl) falls perfectly into the same region of 
GcriT, n) values as a crossover of the temperature depen- 
dence of the diffusion coefficients from weakly dependent 
on T to quite a sensitive function of the temperature. 

A special attention is paid to evaluation of the off- 
diagonal (relative to the site labels) distribution func- 
tions, which describe how fast the loss of the coherence 
occurs. While the study of the generalized diffusion coef- 
ficients (being related to the "velocity- velocity" autocor- 
relation function) allows one to make a conclusion about 
the recrossing phenomenon, behavior of the off-diagonal 
distribution functions shows us how a multiple crossing 
of the dividing surfaces (placed at the neighboring ad- 
sorption sites) by a moving adparticle can proceed. 

Our paper is organized in the following way. In Section 
II we start from the unitary transformed Hamiltonian for 
a dissipative two-level system on a new correlated basis 
17, [m, H^. In Section III using the obtained earlier 
261 non-Markovian equations for non-equilibrium distri- 
bution functions, we investigate a long-time asymptotics 
of the kinetic kernels determining a dissipative motion 
of the adsorbate. A particular case of a completely co- 
herent dynamics, which corresponds either to T — or 
zero coupling limits, is considered. The expression for 
the generalized diffusion coefficients is obtained in the 
next Section. In Section V the Markovian approxima- 
tion for these functions is considered, and a temperature 
behavior of the experimentally observed diffusion coeffi- 
cient is studied in a weak-coupling limit. In Section VI a 
thorough analysis of the critical diagrams separating dif- 
ferent kinds of the adparticle motion is performed, and 
the interconnection of the obtained results with those 
of the previous Section is established. In Section VII 
the off-diagonal non-equilibrium distribution function is 
evaluated at different values of the tunnelling constant at 
assumption of a continuous media; the obtained results 
are considered in the context of the multiple crossing, and 
a conclusion about validity of the Markovian approxima- 
tion is inferred. In the last Section we discuss briefly the 
obtained results and draw flnal conclusions. 
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II. UNITARY TRANSFORMED HAMILTONIAN 
OF THE "SUBSTRATE-ADSORBATE" SYSTEM 

To specify all interactions in the "substrate-adsorbate" 
system we choose a Hamiltonian, which allows site-to- 
site tunnelling of the adsorbate, intrasite oscillations of 
the adparticle between the ground and the 1-st excited 
states within the potential well, and interaction of the 
adparticles with the lattice (both by density and oscilla- 
tion modes). We use the same basic Hamiltonian as in 
Ref. 

Usually, in quantum diffusion problems one can con- 
sider the "substrate-adsorbate" coupling to be arbitrary. 
On the other hand, one-particle characteristics of the sys- 
tem dealt with the intersite hopping and the intrasite mo- 
tion are treated as small parameters. In such a case it is 
useful [13, EE mi to start from the unitary transformed 
Hamiltonian on a new correlated basis, which provides a 
better zeroth-order representation: the sequence of uni- 
tary transformations has the effect of changing to a rep- 
resentation in which the adsorbate is localized at the left 
(L-state) or at the right (i?-state) end of an adsorption 
site, and in which there is a correlated displacement of the 
lattice. Thus, a starting point in our study is the trans- 
formed Hamiltonian of the system of adparticles : 



H = H 



intra ^ -^inier 



Hpp + Hb — H' + Hpp + Hb- (2.1 



The term Hintra describes the lattice-modified intrasite 
dynamics of the adparticle: 



Hintra = y^^ ^'^s(?^s - 1) - f -^^saLasfl, + b-C 



Bs = exp 



(2.2) 



(2.3) 



with a lattice induced operator exponent Bs , where deno- 
tation h.c. means Hermitian conjugation. The other de- 
notations in Eqs. p.2|) - (|2.3p are the following: U means 
the on-site Hubbard repulsion; D, stands for the vibra- 
tional frequency between L- and -R-states with the cor- 
responding creation a^^, or annihilation a^i, Qsr 
operators of the adparticle at a given site s (in these 

notations Ug — a\^asL + ilflOsfl)- It has to be stressed 
that Ml coincides with the oscillation energy between the 
ground and the 1-st excited vibrational states at the as- 
sumption that zero of energy lies midway between these 
two levels [13, [13 • The strength Xsq describes coupling 
of phonons with the energy Hb — X]g ^qb\bq to the os- 
cillation modes of the adsorbate, and only an acoustic 
branch LOq of the substrate motion is taken into account. 
The second summand in (12.11) 



H.nter^ Yjch[B^s^al^as.R + BfJi-ali^a^-L) 

{ss') 

+tpriBi'/ almas' L + Sf/aLa.'fl) (2.4) 



is the intersite tunnelling term with end-changing (end- 
preserving) amplitudes tch (tpr), which are the linear 
combinations tah = ^(ii ± t^) of the initial tunnelling 

pr- 

amplitudes (usually, a condition ti » to is valid, thus an 
underbarrier hopping of the particle is determined by the 
tunnelling constant ti between the excited states rather 
than between the ground states). An abbreviation {ss') 
in (|2.4|) denotes a sum over the nearest-neighbor sites. 

The structure of lattice induced operators B in 
Eq. (|2.4p is similar to that of 



= exp[-^(Af ±(+)5f)(6, - b% 



B^. 



exp[-^(Af ±(-)5f )(&,-&t)], (2.5) 



A' 



' _ 7sg - 7s'9 (±)i^ss' _ Xsq ± Xs' 



huj, 



fkjJn 



where the upper sign corresponds to the first superscript, 
and the strength jsq describes coupling of phonons with 
the density mode of the adsorbate. In a ID case both 
7sg and x^n can be expressed explicitly via the lattice pa- 
rameters [141; for a 2D infinite lattice we shall introduce 
the lattice spectral weight functions [13, [3 to describe 
"substrate-adsorbate" interaction. 
The last but one term of the Hamiltonian (12.11) 



{ss') 

+C^^?insL - nsR)ins'L - n^'j?)} (2.6) 

describes the particle-particle lattice induced interaction. 
For the explicit expressions of the strengths C*"'^/, 
i,j — {O, D}, and explanation of their physical meaning 
the reader is referred to Refs. [13, IS, H^]. It has also 
to be stressed that a direct "adsorbate-adsorbate" 
interaction can be introduced at this stage, and it will 
modify the expressions for C*"'^, to a certain extent. 

To describe the "substrate-adsorbate" interaction we 
consider site-independent end-changing spectral weight 
functions: 

'/Lfl(w)=E [(^''9 ^ 1s'q) + {Xsq + Xs'q)fS{uJ - OJq), 



(2., 



Jrl{uj)^'Y^ [{jsq - 7s'q)~{Xsq + Xs'q)fS{uj - OJg), 
9 

and end-preserving ones: 

Jll{uj)^'Y^ [{jsq - 7s'q) + {Xsq ^ Xs'q)fS{uj ~ UJg), 
9 

jRRi^^) =E [(^«9 ~ 7s'g)~(Xs9 " Xs'q)]'^S{uj - UJg). 



(2.9) 
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The function (|2.7|) describes the intrasite dynamics; 
the functions (|2.8p are related to the intersite end- 
changing processes, while p.9|) are dealt with the in- 
tersite end-preserving processes. At low frequencies the 
end-changing spectral weight functions (labelled by the 
subscript c) are approximately given by 



Jc{uj) 



0, uj < ujo, 

,D-2 



ricUj^ a; > wq, 
and the end-preserving (with the subscript p) ones by 



Jp(w) 



0, 
ripU 



LU > LOq 



(2.11) 



be obtained using the equation for the reduced density 
matrix ps{t) [22l|. Considering the first two terms of the 
Hamiltonian (12.11) as a small perturbation, we can con- 
struct a closed systems of kinetic equations up to the sec- 
ond order in H' . These equations turn out to be nonlocal 
in time, so it is convenient to perform a Laplace trans- 
formation for the diagonal fs.s{t) — {^li^siYs ^^'^ 

i=L,B. 



(2-10) off-diagonal fs^s'{t) 



i=L,B. 



one-particle non- 



equilibrium distribution functions. In the single particle 
limit we obtain 126] the following chain of equations for 
the Laplace transforms fs^s{z), fs,s+n{z) (the index s + n 
means the nearest neighboring site with respect to s): 



with 



Vc = lOG, 



r,p = 12.5G, 



(2.12) 



given in units of the dimensionless coupling constant \Vi 



G 



(2.13) 



In Eqs. (PTTU)) - (PTTO)) D labels dimensionality of the lat- 
tice; M denotes the mass of a substrate atom; uJmax 
stands for the Debye frequency, and the coupling strength 
F is expressed via the mean value of the distortion po- 
tential over the localized Wannier states. 

It is seen from (j2.10|) - (j2.11|) that the lattice is allowed 
to possess a nonzero lowest frequency Wq. At first glance, 
the presence of a gap in the spectrum of acoustic phonons 
on the infinite lattice seems to be rather controversial. 
There were several attempts [H, [l3, HI] to associate 
the value uj^^ with a finiteness of the particle life-time 
To ^ ^^ax i'^ ^ quantum well. We are not going to repeat 
in our article the main reasonings of cited papers, noting 
that the obtained results for the surface quantum diffu- 
sion coefficients are quite insensitive to the value of wq, 
provided the temperature is much higher than huQ/ks 

MM- 

Another interesting feature of the low-frequency de- 
pendence of the spectral weight functions (|2.10p - (|2.1ip is 
that they are scaled like in the case of a bulk diffusion of a 
light inclusion, assisted by optical phonons [12, . The 
only difference is the value of the parameter D, which is 
found to be equal to 5. 

Now let us proceed to the investigation of the adsor- 
bate dynamics using the system of quantum kinetic equa- 
tions. 



III. KINETIC EQUATIONS FOR 
ONE-PARTICLE NON-EQUILIBRIUM 
FUNCTIONS OF ADSORBATE 

The system of quantum kinetic equations for one- 
particle non-equilibrium functions of the adsorbate can 



Winter 

(z) 2/,,,(z)-^/, (^) 



'^(/s,s+n(^) — /s+n,s(^)) 
n 

(3.1) 



^ f s,s-\-n (-2^) fs,s-\-n — 0) — j-^inter 



{I 



(3.2) 



- {linter{z)+ltntra{z)) fs,s+7i{z) + ^~lj;^{z) f s+n,s{^) ■ 

A similar chain of equations has been obtained in 
Ref. [3] , but the authors applied the Markovian approx- 
imation for the kinetic kernels and did not study the 
short-time dynamics of the adsorbate. 

Linearity of Eqs. (|3.1l) - (|3.2p is the result of the single- 
particle approximation, and the question about statistics 
of the adparticle loses its significance. Note, that a simi- 
lar linear approximation could be performed also at low 
coverage of the adsorbate. However, at low-to-moderate 
coverage one has to retain all nonlinear terms and at high 
coverage to include non-equilibrium correlation functions 
into the set of dynamical variables of the abbreviated de- 
scription (2^ . 

Some words have to be said about all constituents of 
the kinetic equations (I3.1l) - (|3.2p . The first terms of r.h.s. 
describe a nondissipative coherent motion of the adsor- 
bate with the renormalized tunnelling amplitude 



1 f , Jpil^) 

-- / dw^V^coth 



rjpk bT sinh 



1 



2ripkBT 



huj 

2kBT 

ripknT 



(3.3) 



In fact, tinter corrcsponds to the polaron band narrowing 
due to the "substrate-adsorbate" interaction. Hereafter 
we use dimensionless frequencies in the units of ujmax and 
temperatures in the units of hujmax/kB- 
The kinetic kernel 

l7nter{z) = ^1Ll{z) + '21Lr{z) + 27iiL(z) (3.4) 
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corresponds to the dissipative intersite motion of the 
adsorbate and describes processes, when the adparticle 
performs series of random site-to-site hoppings (with or 
without the change of its quantum states) owing to the 
interaction with a bath. The kinetic kernel 'jintraiz) 
in Eq. p.2p describes a dissipative intrasite dynamics, 
when the adsorbate during its scattering from the lat- 
tice gets enough energy from the bath to be excited from 
the ground state to the upper level within one adsorp- 
tion site (the opposite process of particle de-excitation 
with a phonon emission is also taken into consideration) . 
The rates jintraiz), Winter (z) can be obtained from the 
Laplace transformation of the kinetic kernels 



: LOr, 



X2 74 



X Re{exp[-((y9:j,(0) - •fx{r))] - cxp[-93^(0)]}, (3.5) 



xRe{exp[-((pii(0) + (^LL(r))]-exp[-(pii(0)]}, (3.6) 
where 



V5x(t) = 



coth 



fvjj 



cos(clit) — i sin(a;r) 



(3.7) 



In Table I we present the amplitudes and the spec- 
tral weight functions Jx{io), relevant to the rates p.5|) - 
p.6p appearing in the kinetic equations p.ip - p.2|) . The 



TABLE L Rates ^x, dimensionless amplitudes A^, and spec- 
tral weight functions Jx{i^) along with their low-frequency 
limits PTU)) - PTT|) . 



A. 



Ja:(^) 



"yintra 

-yLR 



^ f '2cUrnax 

max 

-(ti+to)/2hjj 

max 
max 
max 



Jintraiy) — Jc{i^) 
Jlr{uj) — Jc{u}) 
Jrl{i^) = Jc{l^) 

Jll{(^) = Jp{uj) 

JrrIoj) = Jp{L0) 



functions p.7p yield lattice contributions to the kinetic 
kernels, and the zeroth order Bessel function Jo in p.Sp - 
()3.6p accounts for the particle contribution. The latter 
function ensures a convergence of the time integrals of 
(|3.5p - p.6p at an arbitrary value of the coupling constant, 
though from the mathematical point of view it exceeds 
the required accuracy, being higher than of the second 
order in tunnelling amplitudes. This result is known to 
follow from going beyond the limits of the 2-nd Born ap- 
proximation for the kinetic kernels j2^ . and provides a 
relaxation of the kernels to zero when t — 00. How- 
ever, one has to be careful when dealing with short-time 



dynamics of the system, because in spite of the kinetic 
kernels decay the problem of energy conservation in the 
system appears [231 ■ Due to this reason we will omit the 
Bessel functions in the expressions p.5p - (l3.6p . using in- 
stead a concept of the finite life-time tq — ujq^ of the 
adparticle at a given adsorption cite. 

A study of the long-time asymptotics of the kinetic 
kernels allows one to establish a one-to-one correspon- 
dence between the low-frequency behavior of the spec- 
tral weight functions J (to) and a damping of the ker- 
nels p.5p - (l3.6p at long times. These results are summa- 
rized in Table II, where the constants ai, i — {0, . . . , 3}, 
are introduced just to describe a particular time behav- 
ior of j{t) (in general case, these values are defined 
by the system parameters). It is seen from Table II 



TABLE II: Relation between a low-frequency asymptotics of 
the spectral weight functions and a long-time relaxation of 
the kinetic kernels. 



J(a;)~ 




uj",n>3 


exp(— aoT^) 
exp(— air) 

exp(-a3) 



that we pass from a fast relaxation of the kinetic ker- 
nels at J{oj) ~ w", {n — 0,1}, through long tails at 
71 = 2 to the divergent transport coefficients (which 
are defined via time integrals of p. 5p - (13.61) ) at n = 3. 
Note, that the system dimensionality, as it is seen from 
Eqs. (|2.10p - (|2.1ip . influences the long-time asymptotics 
of the kinetic kernels above all: in the bulk the only lin- 
ear "substrate-adsorbate" coupling does not ensure the 
finite values of the diffusion coefficients. As it has been 
already mentioned, this divergency can be eliminated by 
introduction of the additional channels of the particle 
scattering (electronic friction [3 ; anharmonic terms 
in "adsorbate-substrate" interaction [l^, [l^ or direct 
"adsorbate-adsorbate" interaction), which changes the 
low-frequency asymptotics of the spectral weight func- 
tions. For instance, taking into account the first two fac- 
tors yields damping of the kinetic kernels as exp(— ait) 

Keeping in mind the data of Table II, it is possible to 
evaluate the long-time asymptotics of the end-changing 
(see first three rows of Table I) and end-preserving (last 
two rows) kinetic rates 'yah (t). Both rates in the strong- 

pr 

coupling limit (G > 0.1) decay as Gaussian functions 



7c (t) ~ exp [-77c I In Wo I {ksTT^ + ir)] 



7p(r) - exp 



(3.8) 



(3.9) 



while at the weak-coupling limit (G < 0.01) the first of 
them still decays as a Gaussian function, but the second 
of them behaves as 
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7p(t) - l/T^^^kuT „ ^2r,,fcBT^ ^ ^ 2.15. (3.10) 

We will return to the expressions p.Sp - p.lOp in Section 
V when evaluating a temperature behavior of the diffu- 
sion coefficients. We also leave aside the general cases 
of arbitrary n (both the sub-ohmic, < n < 1, and 
the super-ohmic, n > 1, regimes) in the low-frequency 
asymptotics of the spectral weight functions, referring 
readers to the review [2^. Instead, we consider a limit- 
ing case of the vanishing coupling constant G — >■ 0, or, 
equivalently, the low-temperature limit T — > 0, when the 
adparticle motion becomes completely coherent. For sim- 
plicity, we limit ourselves to the ID, single-band {VI = 0) 
approximation. In such a case the kinetic equations p.ip - 
p.2p can be rewritten as 

fsAt) = 'jt^ J2 (Us+n{t) - UnAt)) . (3-11) 




1 2 3 4 5 6 



FIG. 1: The time dependence of the one-particle nonequi- 
librium distribution functions fo,o{T*) (solid line), /i,i(r*) 
(dotted line), and imaginary part of /o,i(t*) (dashed line), 
given by Eqs. H3.13|) - (|3.14p . The symbol t* denotes a time in 
the units h/ {2to). 



fs,s+n{t)~— —to {fs+n,s+nit)^fs^sit)+fs-n,s+n{t)—fs,s+2n{t)) 

(3.12) 

It is seen that in the site representation we face with a 
coupled chain of equations involving all the lattice labels 
s. However, this system of equations can be solved in 
the wave- vector representation (the details of the solution 
can be found in Appendix), giving the final result via the 
s-th order Bessel functions J, as 



(3.13) 



/,,,+i(t) = Re [fs.s+im + tJs{2toh-^T)J,+i{2tof{2}hl) 

The real part of the off-diagonal distribution function 
can be evaluated by the methods of equilibrium statis- 
tical mechanics. It defines a strength of the transition, 
and does not evolve in time. Instead, the imaginary part 
of the /s_s+„(t) is shifted by the quarter-period with re- 
spect to the diagonal distribution functions (see Fig. 1). 
A particle, initially located, say, at s = 0, starts its mo- 
tion toward the nearest adsorption site s = I. At that 
time the probability to find the adparticle at the site 
s — reduces, the probability to meet the adparticle at 
the site s — 1 increases. The imaginary part of the tran- 
sition probability fs,s+iiT) reaches its maximum value 
when the increase/decrease of the diagonal distribution 
functions becomes maximal (r* w 1). When the parti- 
cle reaches the site s = 1 at r* —2, and the value of 
/i.i(r) is maximal, the inverse motion of the part of the 
wave packet towards the site s — begins, giving the 
negative branch of Im[/o,i(T)] at the subsequent period 
of time. Another part of the wave packet keeps on a mo- 
tion towards the site s = 2 yeilding a positive branch of 
R-e[/o,2(T)] and Re[/i,2(''')] (not presented in Fig. 1). 

In contrast to the classical picture, when a free parti- 
cle performs quasi-continuous motion over the barriers. 



the reason for oscillations of the distributions functions 
is purely quantum mechanical: a superposition of the 
wave packets reflection from the potential barriers and 
tunnelling through them. 

To study the incoherent motion of adsorbate we have 
to obtain expressions for the generalized diffusion coeffi- 
cients. This task is accomplished in the next Section. 



IV. GENERALIZED DIFFUSION EQUATION 
FOR THE ADPARTICLE AT THE SURFACE 

Our further advance is in the manner of the generalized 
collective modes approach [sO, IsH . To find the general- 
ized diffusion coefficient let us solve Eq. p. 21) with respect 
to the hopping probabilities fs^s+n{z) and insert the ob- 
tained result into Eq. p.ip . After grouping the terms, 
one obtains the following equation: 



z/,,,(z)-/,,,(0) 



2t2 



Z + Jinter{z) + 7intra(z) + ithi^) 



~\~^inter \Z 



)} (Y.fs+n.s+niz)-2fsAz)^ 



(4.1) 



''inter"' 



Z + Winter (z) + lintra{z) + 7ll(^) 
xRe i^fs-n,s+n{z) - /s,s+2n(z) 4" fs,s-2n{z)^ . 

The ratio in the braces describes a coherent contribution 
Dcoh{z) = 2tf„tg^/i"^/7tota;(z) to the generalized diffu- 
sion coefficient and can be interpreted \lK [2^ in terms 
of a simple model of band-type motion limited by scat- 
tering from the lattice at temperatures large relative to 
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the bandwidth: 



(4.2) 



where v — atinter/fi denotes the average velocity of the 
adsorbate, a stands for a substrate interatomic spacing, 
and 

ltotal{z) = linter{z) + %ntra{z) + ithi^) (4-3) 

means the total rate of scattering from the lattice. A 
coherent contribution Dcoh{z) characterizes the way in 
which the dephasing limits the band motion of the 
adatom by destruction of the coherence of the hopping 
probabilities /s,s+„(z). Whereas the eigenstate of a free 
particle on the surface is described by a superimposi- 
tion of localized Wannier states (this limiting case cor- 
responds to the ballistic regime of motion) , the coupling 
with the thermal bath induces random fluctuations of 
each phase which destroys the coherence of the state. 

The second term Din(z) = ^mteriz) in the braces of 
Eq. (|4.ip is an incoherent contribution to the generalized 
diffusion coefficient. This is the result expected from the 
random walk model for diffusion with site-to-site hopping 
rate Winter (z), describing processes of the surface phonon 
creation/annihilation when the particle performs a tran- 
sition from one Wannier state to another. 

The last term of the r.h.s of Eq. (|4.ip involves the tran- 
sition probabilities of the particle to perform a long hop- 
ping with |s — s'l > a. In the case of consideration, when 
ti <^ hfl, this term is very small, being of the 4-th order 
in tunnelling amplitude ti. If ti ^ hft, it is comparable 
with the second term of Eq. (|4.ip . being ~ tf. Whatever 
the case, as it is shown in Appendix, the last term does 
not contribute to the overall diffusion coefficient. 

Let us remind that in Section III we assumed the ki- 
netic kernels to be independent of the site label s. This 
assumption leads to the absence of spatial non-locality 
in the expressions for generalized diffusion coefficients, 
so the memory effects only are taken into consideration. 
A spatial inhomogeneity is the subject of separate studies 
(see, e.g. Ref. ^32]). The most general case of the wave- 
vector dependent diffusion coefficients also with the time 
non-locality, that allows one to get a deeper insight into 
dynamics of the system at various time-spatial scales, is 
a challenging topic of non-equilibrium surface diffusion 
theory but lies beyond the scope of the present paper. 

A multiplier ^ fs+n,s+n{z) — 2fs^s{z) at the braces in 

n 

Eq. (j4.ip in the continuous media limit, when the in- 
teratomic spacing tends to zero, converts to the second 
derivative with respect to the space variable (for ID lat- 
tice) times or to the Laplace operator (for 2D lattice in 
absence of the next-to-nearest-neighbor hopping) times 
4a^. It was shown in Ref. [26] that in the continuous 
media limit, and when coupling between the adparticle 
and the surface is strong enough, it is possible to ob- 
tain the Telegrapher's equation for the nonequilibrium 
distribution function n{r,t). This equation is known to 
describe a correlated random walk [33| and, usually, is 



obtained phenomenologically by introducing special re- 
laxation flux terms to the original diffusion equation. In 
(26j it was obtained rigorously by the Markovian approx- 
imation 



DZh{z) 



1 



Z + ltotal{z = 0) 



4 



1 



z + Ifl-^jciz = 0) 



(4.4) 



for the coherent part of the generalized diffusion coeffi- 
cient and a zero-width approximation for the Gaussian 
functions p.8|) - (|3.9|) (which define the incoherent term) 
as a consequence of different time scales for both mecha- 
nisms of dissipation. 

Another remarkable feature of the short-time dynam- 
ics of the adsorbate, described by Eq. (14.11) along with 
approximation (j4.4p for the coherent contribution to the 
generalized diffusion coefficient, follows from the expres- 
sion for a mean square displacement of the adparticle 



lim 



e — ioo 



Dcoh{z) + Dtn{z) 



exp{zt)dz (4.5) 



OJ.2 t-2 

= l^nter{0)t + [exp(-7totaZ (0)i) - 1 + 7tota;(0)i] 

While in the hydrodynamic limit i — > oo Eq. (|4.5p repro- 
duces the Einstein's law for the mean square displace- 
ment 



(Ar(i)2 



= 7mter(0) -|- — 



Qj-2 fc-2 
"''inter"' 



ItotaliO) 



t, (4.1 



in the short-time limit i <C 1 j^totai (0) the second term of 
(|4.5p converts to 4tf^^^j.h~'^t^ . The ballistic term, being 
obtained from the diffusion equation for fs,sit), is quite 
uncommon because in the general case it appears only 
if one uses a Fokker-Planck equation for the distribution 
function, depending on both a coordinate and a velocity 
of the particle [2^, ■ 

At the end of this Section we would like to point 
your attention to the following matter. We call time- 
dependent diffusion coefficients D{t) the generalized dif- 
fusion ones, even though this denotation is usually 
[31I, [13] attributed to the Laplace-transforms of D{t). It 
should be stressed that the time dependence of the ki- 
netic kernels is much more informative than frequency 
one: the generalized diffusion coefficient D{t) is directly 
related to the velocity autocorrelation function C^f (t), 
determined on the adsorption site s. The investigation of 
its temporal behavior can help to visualize the adparticle 
motion both at short and long times, and it is a subject of 
Section VI. As for the next Section, we are going to study 
the temperature dependence of the diffusion coefficients, 
applying the Markovian approximation to Dcoh{z) and 

Din{z). 
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V. TEMPERATURE BEHAVIOR OF THE 
DIFFUSION COEFFICIENTS 

It is known that the experimentally measured values 
of the diffusion coefficients D^xp are usually associated 
with zeroth moments of the generalized diffusion coeffi- 
cients Dq — D{t)dt, which is nothing but the Marko- 

vian approximation Dq = D{z = 0) for their Laplace 
transforms. In fact, in experimental conditions one deals 
with evaluation of the mean square displacement of the 
adparticle at times much larger than l/^totai{0)- Thus, 
measuring (Ar(t)^) one has to be sure that the influence 
of the transient states is excluded, and duration of the 
atom-tracking procedure is large enough to fall into the 
hydrodynamic region t ~> od. Otherwise, the value -De^p 
will differ from its theoretical prediction We con- 
sider the transition regimes of the adparticle motion in 
the next Section. In this Section we investigate the tem- 
perature behavior of the experimentally observed diffu- 
sion coefficients in the framework of two-level dissipative 
model. 

First of all, let us determine the conditions of the va- 
lidity of Markovian approximation. It is believed @, 
that memory effects can be neglected if the time scales 
describing the adsorbate motion and those of the lattice 
dynamics are well separated, ujmax/^ ^ 1- Stronger 
substrate-adsorbate coupling favors the Markovian ap- 
proximation, while the weak-coupling limit usually re- 
quires consideration of the memory effects at the initial 
stage of the adparticle motion. 

If the Markovian approximation is valid, the diffusion 
coefficient is determined by the Einstein's law (|4.6|) as 



7mter(0) 



oj-2 t.--^ 
"''inter"' 

ItotaliO) 



(5.1) 



Taking into account the expressions p.3p for the renor- 
malized tunnelhng amplitude and Eqs. p.5p - p.6l) . (|4.3p 
for the kinetic kernels, and noting their time dependence 
(I3.8p - p.f0p . it is easy to perform an integration over r 
and to obtain the final result. 

Thus, in a strong-coupling limit G > O.l the diffusion 
coefficient is completely defined by the incoherent term: 



^strong — 



.(0) 



which, in its turn, is being determined by the end- 
preserving processes. 

In the weak-coupling limit G < O.OI, and at a reason- 
able assumption ti <^ hil, the diffusion coefficient 



D 



weak 



2 r 



8exp(-y)p(0)) j 
(7p+7c) + 15.3) 



contributions. The end-changing 7c and the end-preserv- 
ing 7p rates can be presented in the following form: 




is determined by both incoherent (the first term of 
Eq. (|5.3p ) and coherent (the second term of Eq. (|5.3p ) 



7c 



TT exp(-77c| lna;o|/4fcBT) 



7p 



^0 



VkBT 

(A-2)ri^kBT 
^ 

- ArjpkBT 



(5.4) 



(5.5) 



It has to be noted that the expressions for 7c in both 
limits coincide with those of Ref. 18], while the expres- 
sion for 7p in the weak-coupling limit differs from the 
result of cited paper, which was obtained as the mul- 
tiphonon expansion of the corresponding end-preserving 
rates. Keeping in mind the power law behavior p.lO|) 
of the end-preserving kernel at weak-coupling, one can 
justify that a multiphonon expansion is not valid in this 
particular case. 

Now we have all necessary conditions to evaluate the 
temperature behavior of the surface diffusion coefficients. 
In Fig. 2 we present the log-log plot of the diffusion co- 
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FIG. 2: Temperature dependence of the diffusion coefficient 
in a weak-coupling limit. Model parameters: G — fO~'^, 
LJo /i^max — 5 X fO"*, ti / hijjmax = fO~^. Lines with triangles, 
diamonds, stars, circles and squares denote, correspondingly, 
the values of Q./u)max=Q-2, 0.25, 0.3, 0.5 and f . 

efficients in the weak-coupling limit as functions of tem- 
perature. The sum of Dcoh and Din gives rise to quite 
a complex temperature behavior of the overall diffusion 
coefficient, but for all values of VL shown, D^eak is a rel- 
atively insensitive function of the temperature. 

The arguments about validity of the Markovian ap- 
proximation at values of fi, presented in Fig. 2, are not 
contradicting with the condition ojmax/^ ^ 1, as long 
as we are not interested in the adsorbate motion at short 
time scales: according to Eq. (|4.4p . the zeroth moment of 
Dcoh{t) can be calculated at any value of SI. We will see 
in the next Section that it is not true if one investigates 
the intermediate regimes t ^ ^/^totai{^), and it is neces- 
sary to consider a diffusion equation which is nol-local in 
time. 
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In Fig. 3 we show the plot of the full diffusion coeffi- 
cient at stronger coupling to the lattice but still in the 
weak-coupling regime. 




FIG. 3: Temperature dependence of the diffusion coefficient 
in a weak-coupling limit. Model parameters: G = fO~^, 
LOo/ujmax = 5 X 10"^, t\/huj,-nax ~ fO~^. Lincs with trian- 
gles, squares, circles and stars denote, correspondingly, the 
values of n/aJmaa;=0.25, 0.5, f , and 2. 



The temperature dependence is quite different from that 
shown in Fig. 2: at low temperatures the diffusion is dom- 
inated by the coherent contribution, which is a strongly 
decreasing function of T. According to (j5.3|) - (l5.4l) . the 
coherent contribution to the diffusion coefficient behaves 
as exp^Ec/ksT) at low temperatures and as \/kBT at 
higher temperatures. At high temperatures the main 
contribution comes from the incoherent term, and D^^ak 
is a slowly increasing function of T. 

Let us briefly discuss the temperature behavior of 
Dweak at low tcmpcratures, shown in Fig. 3. A qual- 
itatively similar increase of the diffusion rate of H on 
Cu(OOf) below 20 K was observed experimentally in 
Ref. However, in our case this behavior is just a 
result of the used two-level dissipative model, whereas in 
Ref. 3 it was attributed to the change of nonadiabatic 
response of the thermally excited electron-hole pairs to 
the diffusing particle. When plotting Figs. 2 and 3 we just 
emphasize that a crossover from one kind of temperature 
dependence of the diffusion coefficient to another takes 
place even in a relatively simple model. In the next Sec- 
tion we will show that this crossover coincides perfectly 
(regarding to the coupling constant G) with the change 
of the character of adparticle dynamics at short times, 
when the memory effects have to be taken into account. 



VI. TRANSITION REGIMES OF THE 
GENERALIZED DIFFUSION COEFFICIENTS 

We have already mentioned that for visualization of the 
processes of adparticle motion at short and intermediate 
times it is much more convenient to perform an inverse 
Laplace transform of the generalized diffusion coefficients 
according to 



Dcoh{t) =Re 



aty 



f 1 

lim / dz exp{zt) 

<^-^0 J Z + ^total 



(^) 



Re 



OO ^ 

^exp(-z,t)— — 

.i=l ^ 'to 



3.1) 



The summation in (j6.ip in accordance with the residue 
theorem runs over all poles Zi of the integrand, which 
obey the condition Re[zi] < 0. The expression (|6.II) re- 
sembles the results of the generalized collective modes 
theory [s^, HH , postulating the additive contribution of 
each collective excitation to a particular time correlation 
function. In our case, the summation is extended to the 
infinite number of poles, and a main contribution comes 
from the terms with maximal values of Re [2:1] and weight 
factors + %otai('^i)]~'^ ■ ^^^'^ mentioned that 

the expression for Dcoh (t) can be even more complicated 
if one deals with poles of the order n > 1. 

The expression for Din (t) follows immediately from the 
definition (j3.4p of the kernel jinteriz) and can be writ- 
ten down via the end-changing/end-preserving functions 
(13:81)- (I3l0ll as 



= {ah/h)^Re [7,(t) + 7p(i)] 



(6.2) 



We evaluate the generalized diffusion coefficients in the 
weak-coupling regime when the memory effects are im- 
portant at the initial stage of the adparticle motion. 
In Fig. 4 we present the time dependence of D{t) — 
Dinit) + Dcoh{t) at different temperatures. One can draw 
some conclusions after observation of the plots. 

First of all, we see a non-monotonic behavior of D{t), 
and these oscillations become more pronounced when 
the system temperature decreases. At low temperatures 
thermal fluctuations of the lattice are very small. There- 
fore, the lattice distortion caused by interaction of the 
adsorbate with the phonon subsystem has no time to 
relax, and the initial profile of the lattice potential has 
no time to be restored after the particle passage. As a 
result, the effective barrier is higher than its adiabatic 
value, and the particle starts to oscillate being caged in 
the deformed potential well. Such a behavior of D{t) 
is observed even at the temperatures comparable with 

Secondly, at temperatures ksT /Tiujmax < 0.5 (which 
are not presented in Fig. 4) these oscillations persist on 
the time scales, which are by two orders of magnitude 
higher than the inverse Debye frequency. But even at 
higher temperatures there are evident memory effects. 
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FIG. 4: Generalized diffusion coefficients calculated at G = 
10""*, t\/fkjjmax ~ 10~^, Q./ijJmax = 1, aud different tem- 
peratures: kBT/hijJmax=0.5 (triangles), 0.75 (circles) and 1 
(squares). 



In Ref. [26| we calculated the generalized diffusion co- 
efficients at different values of T and ft. The general 
tendencies observed in [2^ can be formulated as follows: 
a low temperature T and a high vibrational frequency 
favor the oscillation dynamics of the particle, and so does 
a weak "substrate-adsorbate" interaction. 

It would be interesting to relate this non-monotonic 
behavior of D{t) to a possible recrossing phenomenon 
0, B HEl- One can attribute the negative branches of 
D{t) to the backward motion of the adsorbate: the par- 
ticle may cross the dividing surface, located at the ad- 
sorption site s, due to the lattice distortion that "pushes" 
the particle in the opposite direction (with respect to that 
of the initial instant of motion). The nature of oscilla- 
tions of D{t) is different from that of fs,s{t), presented in 
Fig. 1. While in the coherent regime the only reason of 
the non-monotonic behavior is an interplay between the 
processes of transition and reflection of the wave packet, 
associated with the adparticle, the non-monotonic inco- 
herent motion is determined by the adsorbate scattering 
on the substrate atoms (one can verify that only Dcoh {t) 
contributes to the oscillatory adpartacle dynamics). 

Now let us ask the question: what happens when one 
increases the value of coupling constant G? An intu- 
itive answer would state that oscillations of D{t) disap- 
pear at a moderate-to-strong coupling. Indeed, at strong 
coupling, when the energy exchange between the parti- 
cle and the substrate atoms is faster, one can apply the 
Markovian approximation ()4.4p for the coherent contri- 
bution to the generalized diffusion coefficient, which leads 
to the exponential relaxation of Dcohit). However, fine 
features of such transition regimes, when the character of 
the adparticle motion changes from oscillatory to mono- 
tonic, need a thorough analysis in the framework of the 
non-Markovian approach. 



Thus, if one increases the coupling constant G un- 
til the oscillations of D{t) disappear, one can obtain a 
critical value Gcr{T, ft) as a function of the temperature 
and vibrational frequency, which separates two regimes 
of time evolution: there is a plain relaxation of D{t) 
at G > GcriT,n), and a non-monotonic behavior at 
G < GcriT^ft). A detailed analysis shows that, at first, 
the negative branch of D{t) rises over the time axis, at 
that oscillations of the generalized diffusion coefficients 
still remain. So recrossing is vanishing, but the particle 
moves toward the nearest adsorption site at one moment 
faster, at another slower as if meeting obstacles. 

At higher values of G the oscillations completely dis- 
appear, and the adsorbate motion is governed mainly 
by the incoherent term Din{t). At G ^ 0.1 no coher- 
ent contribution is evident. Moreover, in the expres- 
sion (|6.2p for Din{t) one has to use the strong-coupling 
limit for the end-preserving kernel p.9p rather than the 
weak-coupling form p.l0|) . At such values of G the end- 
preserving rate jp{t) solely defines the adparticle dynam- 
ics, and the temperature dependence of the diffusion co- 
efficient is given by Eq. (|5.2p . 
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FIG. 5: The critical values of the coupling constants 
Gcr{T,Q*) at ^1* /ujmax ~ 1, as functions of temperature T, 
that separate domains of the monotonic adparticle dynamics 
(A), oscillations without recrossing (B), and an eventual re- 
crossing (C). The inset shows the temperature dependence of 
Gcr{T,Q*) that separates the strong-coupling (1) and weak- 
coupling (2) domains. 

We present the above mentioned transition regimes of 
the adparticle motion in Fig. 5 and 6 as plots of Gcr 
vs. temperature (at the fixed value of vibrational energy 
n* /iOmax = Ij Fig. 5), and vs. vibrational energy (at the 
fixed value of temperature ksT /humax = 0.5, Fig. 6). 
When inspecting these curves, one can observe a remark- 
able feature: the transition domains of different dynamic 
regimes of D{t) with respect to the coupling constant 
G coincide with the regions, where the temperature be- 
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FIG. 6: The critical values of the coupling constants 
Gcr{T* ,Q) at ksT* /hwmax ~ 0.5, as functions of the vibra- 
tional frequency Q., that separate domains of the monotonic 
adparticle dynamics (A), oscillations without recrossing (B), 
and an eventual recrossing (C). The inset shows the frequency 
dependence of Gcr{T* ,0) that separates domain (1), where 
Dcohit) decays as an exponential function (see Eq. (|4.4|l '). from 
that with essentially non-Markovian behavior of Dcohit). 



havior of the diffusion coefficients changes from a weakly 
dependent function of T to quite a sensitive function of 
temperature. Namely, ttie region witfi a recrossing (the 
C-domain in Fig. 5) maps to a weakly dependent temper- 
ature regime (see Fig. 2), and the region with monotonic 
motion of the adparticle (the A-domain in Fig. 5) corre- 
sponds to a strongly dependent temperature regime (see 
Fig. 3). The domain B in Fig. 5, evidently, corresponds to 
transition of the temperature behavior of experimentally 
measured diffusion coefficient from weakly to strongly de- 
pendent function of T. 

The same tendencies can be traced in Fig. 6, where 
the plot of Ccr vs. Vl is presented at the fixed temper- 
ature. High values of the vibrational frequency extend 
the domains (B and C) of the non-monotonic adparti- 
cle dynamics. The insert in Fig. 6 shows that a tran- 
sition from essentially non-Markovian dynamics to the 
case, when Dcohit) decays as an exponential function (see 
Eq. (j4.4[) ). occurs in a moderate-to-strong coupling do- 
main. Of course, the most informative would be a three- 
dimensional plot GcriT, fi), but even such sections of the 
critical coupling constants at fixed temperature and vi- 
brational energy, as in the above presented figures, give 
much food for thought. 

One of the assumptions inferred can be formulated 
as follows. If the tendencies, presented in the "critical 
diagrams" in Figs. 5 and 6, remain in more sophisti- 
cated models describing quantum surface diffusion, and a 
correspondence between the "recrossing-monotonic mo- 
tion" transition, and the change in the character of a 



T-dependence of the experimentally measured diffusion 
coefficients D^xp is valid in general, then we could give 
a prognosis about the temperature behavior of D^xp, lia- 
ving only the information about the dynamics of the ad- 
sorbate at the initial stage of its motion. Otherwise, 
one has to measure the mean square displacement of the 
particle at times large enough to be sure that all rem- 
nants of the transition regimes are excluded. As it has 
been already mentioned, due to comparatively slow pro- 
cesses of quantum diffusion at low temperatures (where 
the aforesaid transition regimes are the most pronounced 
and durable), it would be useful to reduce the time of 
atom-tracking experiment, because the longer action of 
the measuring instrument, the greater influence on the 
system occurs. 

It should be also noted that investigation of the tran- 
sition regimes in the framework of two-level dissipa- 
tive models is a topical problem [2^, and much efforts 
are put into a study of the system crossover from one 
kind of its dynamics to another. In a recent paper [35 1 
the authors obtained the "coherence-incoherence" tran- 
sition diagrams for the dissipative two-level system with 
a nonzero bias and a sub-Ohmic bath as functions of the 
power index n in a low-frequency asymptotics J(w) ^ uj" 
of the spectral weight functions. Though their results 
were obtained in zero-temperature limit, it would be in- 
teresting to generalize such a model for T ^ case to 
investigate whether there is any other reason (except the 
G ^ limit) for the "coherence-decoherence" transition. 

At the end of the present Section we would like to 
touch upon the study of the recrossing phenomenon once 
more. Though the considered two-level dissipative model 
is too simplified to describe real "substrate-adsorbate" 
dynamics, especially in a comparison with the models 
traditionally used during a direct evaluation of the "flux- 
flux" time correlation functions [2^, it is worthy to 
compare them. On the one hand, in our model the 
lower temperature, the more pronounced oscillations of 
the generalized diffusion coefficient observed. This is dif- 
ferent from the results of Ref. Q where high T favors 
the recrossing phenomenon. On the other hand, if we 
are in the low-temperature weak-coupling domain and 
increase the value of G, we will observe that there is a 
gradual transition from the domain C with a pronounced 
recrossing to the domain B (see Figs. 5 and 6) with os- 
cillations but without recrossing. On the contrary, in the 
high-temperature region with ksT /fujjmax ~ 1 a gradual 
increase of the coupling constant will give us a sudden 
transition of the coherent part Dcoh{t) of the general- 
ized diffusion coefficient from the domain C to the re- 
gion A of the monotonic adparticle dynamics. However, 
this feature (not presented graphically because we inves- 
tigate the crossover for the overall generalized diffusion 
coefficient D{t) rather than for Dcoh{t)) is smeared out 
by the contribution of the incoherent term, which has a 
"long tail" behavior (|3.10p . If one could eliminate these 
long tails, introducing additional channels of the adpar- 
ticle scattering, one would obtain an interesting result: 
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though the recrossing is less pronounced at high T, it per- 
sists in a wider domain of couphng constant (the region 
B completely vanishes). Such a behavior would resemble 
the results of Ref. [s']. It has to be noted that the pre- 
sented model in the general features corresponds to the 
model F of the cited paper (no surface motion is possi- 
ble, both thermal fluctuations and lattice distortions are 
permitted) . 



VII. TIME DEPENDENCE OF THE 
TRANSITION PROBABILITIES fs,s+n{t) 

In this section we study the off-diagonal non- 
equilibrium distribution functions fs^s+n{t)- There are 
some reasons to look closer at the short-time dynamics 
of the above mentioned transition probabilities. First of 
all, the rate of the decay of fs,s'{t) defines the time scales 
at which the contribution of Dcoh (t) vanishes (let us re- 
mind that to derive the expression for the coherent term 
of the generalized diffusion coefficients we solved Eq. p.2p 
with respect to fs^s'{t), and inserted the obtained result 
in the upper equation p.!!) ). On the other hand, the off- 
diagonal non-equilibrium distribution functions fs^s+n{t) 
can be related to the time dependent "flux-flux" cross- 
correlation functions, determined at the adjacent sites s 
and s + n. These time correlation functions are known to 
describe a multiple crossing or a multi-hopping regime 
0, Q . The multi-hopping facilitates an increase of the 
transport coefficients, in contrast to the recrossing, which 
reduces the total rate of the adparticle escape and lowers 
the value of the diffusion coefficient. 

We evaluate the non-equilibrium transition probabili- 
ties fs^s+n{t), solving Eq. p. II) with respect to the diag- 
onal distribution functions fs.sit) and inserting the ob- 
tained result in Eq. p.2p . The final expressions for the 
real and imaginary parts of fs.s+n{t) are 

Re[/, s+n{z)] l£:£+lM +o> 

Z + Re [jinteriz) + Jintra{z) - JLhi^)] 



of (|7.2I) in the limit G — ^ gives us the expression p.l4p 
for lm[fs^s+n{t)]. 

To study the time dependence of the transition prob- 
abilities we perform a certain simplification of Eq. (|7.2p . 
Namely, we pass from the difference form in the numera- 
tor of (|7.2p to the derivative with respect to the adatom 
coordinate r. In fact, we perform the continuous me- 
dia approximation, similarly like it was done in Ref. [26} , 
when we obtained a Telegrapher's equation, describing 
the correlated random motion of the adparticle. After 
such a transformation, we are not bounded to a spe- 
cial geometry of the lattice any more and can rewrite 
Eq. (17. 2[) in the wave- vector-frequency representation for 
the imaginary part of the transition probability P{k,z) 
as follows: 



Im[P(fc,z)] 



hk 


Dcoh{z) 






Dcoh{z) + Din{z) 



(7.3) 



Now, performing inverse Fourier and Laplace transforma- 
tions, we can evaluate the (r, t)-dependence of the imag- 
inary part of the transition probability. In Figs. 7 and 8 
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FIG. 7: Time dependence of the normalized function 

Im[P{r*,t)]/[lm[P{r',t)]Ua. at r* = r/a = 1, G = 10~\ 
fl/ujmax = 1, ksT/hjJmax ~ 1, and tunnelling amplitude 
— 10~^. Solid and dashed lines correspond to the 



where only linear terms in the tunnelling amplitude ti are 
retained. It has to be stressed that in the single particle 
limit there is no coupling between the real and imaginary 
parts, while in the general case of nonlinear kinetic equa- 
tions Re[/s,s+„(2:)] and lm[fs^s+n{z)] are coupled to each 
other. The second conclusion, which is to the point, is the 
following: the decay rate of Re[fs.s+n{z)] is almost the 
same as that of the coherent part Dcoh{t) of the gener- 
alized diffusion coefficient, which is defined by the kernel 

Itotal {z) = ^inter{z) + ^intra{z) + itb^z)). ObvioUSly, 

in zero-coupling limit Eq. (j7.1|) reproduces the value of 
fs,s+n{t) which does not depend on time in the coherent 
regime of motion, and the inverse Laplace transformation 



Markovian and non-Markovian approximations. 

we present the time dependence of the above mentioned 
functions (normalized at their maximum values) at r* — 
rja = 1. The evaluation is performed for two different 
values of the tunnelling amplitude ti/hwmax = 10~^ and 
10^^. In the first case we observe that the adparticle 
reaches the point r* at times of about 10'' inverse De- 
bye frequency. The adparticle motion is very slow, the 
lattice has enough time to relax at such huge times, and 
we do not observe any noticeable difference between the 
Markovian approximation for transition probability and 
the case, when the memory effects are taken into account. 
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FIG. 8: Time dependence of the normalized function 

Im[P{r*,t)]/[Im[P{r*,t)]Ua. at r* = r/a = 1, G ^ 10-^ 
0,/uJmax ~ 1, ksT /lnomax = 1, and tunnelMng amplitude 
ti/tujmax = 10 Solid and dashed lines correspond to the 
Markovian and non-Markovian approximations. 



A situation changes drastically (Fig. 8), when the value 
of the tunnelling amplitude increases by three orders of 
magnitude. The decay time of P{r*,t) is then compara- 
ble with that of the generalized diffusion coefficient (see 
Fig. 4 for comparison). The adparticle arrives at the 
point r* at times of about seven Debye periods 2n/uJmax-, 
when the lattice has not relaxed completely, and the 
memory effects arc still important. Thus, a transition 
state that originally obstructs a multi-hop from the site 
s to the nearest neighboring one, is stabilized by the lat- 
tice fluctuations: the height of the effective barrier at 
s + n can be lower than its adiabatic value at some in- 
stant of time, allowing the particle to perform a multiple 
crossing. We would like to note the strongly aperiodic 
behavior of the transition probability in contrast to the 
smooth relaxation in the Markovian case. The negative 
branches of P{r* , t) mean that the "transition window" 
for the multiple crossing is not always open, and at other 
instants of time the effective barrier height can exceed its 
static value, prohibiting multi-hops of the particle. On 
the other hand, the width of oscillations growths in time 
while their amplitude reduces. It means that multiple 
hops are less probable when time increases, but the pe- 
riod favorable for them to proceed becomes longer with 
increasing t. 

We do not present plots for the real part of transi- 
tion probability, because its time behavior does not differ 
from that of D{t) (see Fig. 4). At the same time we have 
to emphasize that all the curves in Figs. 7 and 8 corre- 
spond to the normalized values of transition probabilities, 
while the non-normalized values in the continuous media 
limit would be much smaller than those of the distribu- 
tion function n{r* , t) at the starting point r* = of the 
adparticle motion. In fact, n(r* , t) tends to infinity at 



r* = 0, but even in the coherent limit and for a discrete 
lattice model the curve for lm[/o,i(i)] lies much lower 
(see Fig. 1) than that for /o,o(0- I* is obvious that the 
higher transition probability, the greater contribution of 
the multi-hops to the overall diffusion coefficient. 

Thus, only rigorous evaluation of the site-dependent 
transition probabilities fs,s+n{t) for a lattice with given 
geometry would remove possible contradictions between 
the results, obtained for the real discrete structure of the 
surface, and our model case of a continuous media. Nev- 
ertheless, we believe that the basic conclusions drawn in 
this Section arc quite reliable in the context of a qualita- 
tive analysis of the adparticle motion. 



VIII. CONCLUSIONS 

In this paper we made the systematic analysis of the 
dynamics of the adparticle, which performs an under- 
barrier tunnelling from one adsorption site to another 
as well as the vibrational transitions between two levels 
inside a quantum well, when its motion changes from co- 
herent to incoherent due to interaction with acoustic sur- 
face phonons. We investigated the short-time dynamics 
of the particle in the case when its vibrational frequency 
O is comparable with the Debye frequency uJmax, and 
memory effects have to be taken into account. Though 
the latter are traditionally thought to be important only 
at transition regimes, when the adparticle motion is nei- 
ther ballistic nor purely diffusive, and the experimentally 
measured diffusion coefficients do not have to remember 
the details of intermediate period, we showed that there is 
a close relation between the adparticle dynamics at times 
r lOw" J,^ and the temperature dependence of diffusion 
coefficients. Namely, as the coupling constant increases, 
the adparticle motion (initially oscillatory) becomes more 
and more smooth indicating that the temperature behav- 
ior of the diffusion coefficients D{T) should change from 
weakly dependent on T to quite a sensitive function of 
the temperature. Though we restricted ourselves by a 
comparatively simple two-level dissipative model, we be- 
lieve that the above mentioned correspondence between 
the short-time dynamics of the adsorbatc and the tem- 
perature behavior of the diffusion coefficients is valid for 
more sophisticated systems. If so, this "T vs. f cor- 
respondence could be helpful at experimental evaluation 
of the diffusion coefficients because it would allow us to 
give a prognosis about their temperature behavior, ha- 
ving only an information about the adatom dynamics at 
the initial stage of its motion. 

We showed that the coherent term of the generalized 
(time dependent) diffusion coefficient, which is defined by 
the adparticle scattering from the lattice, is responsible 
for the rccrossing phenomenon at weak-coupling, but at 
high-coupling regime its contribution diminishes, and the 
particle motion is completely determined by the incoher- 
ent term. We performed a quantitative analysis of such 
transition regimes in terms of the critical coupling con- 
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stant GcriT,n) which depends on the temperature and 
vibrational frequency. We compared the temperature de- 
pendence of the recrossing with the results obtained by 
a direct calculation of the "velocity-velocity" quantum 
time correlation functions 0, @ , and analyzed both com- 
mon and different features of the adparticle dynamics. 

While the generalized diffusion coefficients are con- 
nected to the "velocity-velocity" autocorrelation func- 
tions, the transition probabilities (non-diagonal distri- 
bution functions fs,s'{t)) can be related to the cross- 
correlation function "velocity-velocity" . A study of the 
transition probabilities is two-fold interesting: on the one 
hand, it shows how fast a loss of the adparticle coherence 
occurs; on the other hand, it allows to draw a conclusion 
about the phenomenon of multiple crossing. The con- 
tribution of multiple jumps to the diffusion coefficient is 
enhanced by the fact that in a double jump the random 
walker goes a double distance. Since the diffusion coeffi- 
cients is defined by the squared jump length, the double 
jump can be important even if its jump rate is compara- 
tively smaller than that of a single jump. 

We investigated the time dependence of the transition 
probabilities in the model case of continuous media. It 
was shown that at very small values ti/Hujjnax = 10~^ 
of the tunnelling amplitude the real part of the tran- 
sition probability (which defines a multi-hop strength) 
decays at the same rate as the generalized diffusion coef- 
ficients, while the particle approaches the nearest adsorp- 
tion site at much later times (which corresponds to the 
maximum value of the imaginary part of fs,s'(t)). The 
particle moves very slow, the lattice has plenty time to 
relax, and the Markovian approximation is quite appli- 
cable if one studies non-equilibrium transition probabil- 
ities. Contrary, at ti/huj„iax = 10^^ the typical times 
of decay for both real and imaginary parts of the tran- 
sition probability are of the same order, and remnants 
of the memory effects have an influence on the multiple 
crossing, which is characterized by an aperiodic oscilla- 
tory function. We also verified that the zero-coupling 
limit reproduces the results for a coherent motion, when 
the strength of transition Re[fs,s+n{t)] does not change 
in time, and lm[fs,s+n{t)] is expressed analytically via 
Bessel functions. 

The memory effects could be said to preserve a parti- 
cle coherence to a certain extent: a portion of relaxation 
energy of the lattice is delivered to the adparticles pre- 
venting them from thermalization and maintaining the 
transition regimes from the coherent (ballistic) motion 
to the incoherent (diffusive) one. This reasoning would 
be even more realistic if one introduced the additional 
channels of the adparticle scattering: an electronic fric- 
tion I1 4| . [ 2011 ■ a non-linear "adsorbate-substrate" inter- 
action or a direct "adsorbate-adsorbate" interaction 
[2l| . First of all, this would solve the "long-tails" problem 
of the kinetic kernels [l^ and ensure their convergence 
without introduction of any additional parameter like a 
particle life-time in a quantum well. On the other 
hand, additional interactions introduce new typical time 



scales, which could be well separated (the Markovian pic- 
ture is then valid) or close to each other (then the non- 
Markovian approach is necessary). Besides, by taking 
an "adsorbate-adsorbate" interaction into account one 
can go beyond the limits of small coverage and study a 
concentration dependence of the diffusion coefficients in 
addition to their temperature behavior [36*1 . We believe 
that all these directions are quite interesting from a view- 
point of the study of transition regimes of the adsorbate 
and could be the subject of future investigations. 
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Appendix 



To solve Eqs. p.lll) - (|3.12p let us perform at first a 
Fourier transformation for the creation/annihilation op- 
erators 



, ex'p(iks)a 



as = — ^ exp{-iks)ak, 
(A.l) 

passing from the site representation aj, to wave- 
vector representation a|,, Uk with k =^ [2t: /N)m [k = 
{2Ti/N){m + 1/2)] for the lattice with the odd [even] 
numbers of adsorption sites and m = —N/2,— N/2 + 
l,...,iV/2-l [if A^is even] or m = -(A^ - l)/2, -(iV - 
l)/2 1, . . . , (A - l)/2 [if N is o dd]. 

Inserting (jA.ll) in Eqs. (|3lT|) - (|XT2)) for the coherent 
motion, we obtain the following equation for the inter- 
mediate distribution function Fk_k' (t) 



{alakY: 



dt h 
which is easily solved, giving 



(cosfc- cosfc')-Ffc,fc'(t), (A.2) 



Fk,k'{t) = exp 



-2ito 



cos k — cos k')t 



(A.3) 



In the infinite lattice limit A^ — > cxd, we can pass from 
summation over k to integration over continuous wave- 
vector according to J2k' ' ' — ^ 5¥ /-tt " 'dk. After 
that we obtain an integral representation for the diago- 
nal one-particle non-equilibrium distribution function as 
follows: 



dk[coa{ks) + i sin(fcs)] 



X exp 



-2ito 



cosk 



j2(2ion-V), (A.4) 
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which is nothing but the squared s-th order Bessel func- 
tion. Similarly, we can obtain the expression p.l4p for 
the imaginary part of the transition probability fs^s+i{t)- 
To show that the last term in Eq. (|4.ip does not con- 
tribute to the mean square displacement (Ar(t)^) — 



necker delta-symbol, which in the infinite lattice limit 
converts to Dirac delta- function S{q), and integrating 
(jA.6|) by parts with taking into account (jA.5| . one can 
verify that only the first term in braces of Eq. (lA.sp con- 
tributes to the evolution equation for {Ar{t)^): 



N 

,2 „2 



J2 s^fs,s(t) of the particle let us perform for simplicity 

the Markovian approximation for the generalized diffu- 
sion equation (I4.ip . Using Fourier transformation (jA.ip 
we obtain the evolution equation for the intermediate 
distribution function in the following form: 



d{Ar{tY 
dt 



= [Dcoh{Q)+D^nmh / dk 



dt a2 



{2[D,oh{0) + A„(0)][1 - cos(fc - k')] 



XI dq[l - COS q]Fk,k-q{t)^S{q) = 2[D,oh{0) + D,i$]J\] 



+ Dcoh (0) [cos 2k + cos 2fc'-2cos(fc + fc')]}, 

where the second term in the braces is related to the last 
term in the r.h.s of Eq. (I4.ip . which involves the "long 
distance" transition probabilities fs±n.s^n{t), fs,s^2n{t)- 
The evolution equation for the mean square displace- 
ment can be written down as follows: 



while the second term vanishes at the integration over k: 
(A.5) 



d{Ar{ty 
Jt 



N 



= ^ ^ exp{iqs)Fk,k-q{t) 



s=l k.q 



= -«'E^^-.^-'^Wd^ \J2eM^1^)j ■ (A.6) 
Noting that the sum in the brackets yields N times Kro- 



J dk J dq {cos{2k) ^ cos{k — q) — 2 cos{2k — q)} 

— TT — TT 

TT 

r 

xFk,k-q{t)-^6{q) = j dk[2 cos{2k) ^ COS k]FkMM 



because Fk^k{t) = 1, and all derivatives of Fk^k{t) with 
respect to wave-vector vanish. 

The presented above calculation can be generalized to 
the 2D case, or to the case when the memory effects are 
taken into account. As for the 2D coherent regime, the 
only modification of Eqs. p.l3|) - p.l4|) consists in dou- 
bling of the power indexes at Bessel functions. 
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